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Abstract 



The presence of resonances modifies the passage of light or of electrons through 
a disordered medium. We generalize random matrix theory to account for 
this effect. Using supersymmetry, we calculate analytically the mean density 
of states, and the effective Lagrangean of the generating functional for the 
two-point function. We show that the diffusion constant scales with the 
effective mean level spacing. The latter exhibits a resonance dip. These 
facts allow us to interpret experimental results on light scattering for different 
concentrations of resonant scatterers. 



I. INTRODUCTION AND MOTIVATION 

During the last decade, many phenomena related to the passage of waves through a 
disordered medium have been intensely studied. This holds both for electrons, i.e. for the 
amplitude waves of the Schrodinger equation U, and for light, i.e. for classical electromag- 
netic waves 0. As examples, we mention conductance fluctuations in mesoscopic probes, 
the analogous effect of speckle patterns in the transmission of light, and weak and strong 
localization effects including the enhanced backscattering of electrons or light from a disor- 
dered medium. For electrons, disorder is typically caused by impurities. For light, disorder 
is often produced artificially: A powder is immersed into some liquid. Random scattering 
occurs if the indices of refraction of powder and liquid are sufficiently different. 

The present paper has been motivated by an experiment performed several years ago 
H on the scattering of light by a disordered medium. In this experiment, the powder used 
(T1O2) consisted of grains with a size distribution centered around a diameter of 220 nm. For 
such granules, a Mie resonance occurs close to the wavelength A ~ 630 nm of the laser light 
used in the experiment. This led to a resonance enhancement of the diffusive scattering. The 
effects of such enhancement were seen by comparing the diffusion constant T> (determined 
from the intensity autocorrelation function versus frequency of the transmitted light), with 
the transport mean free path I (determined either from weak localization, i.e. from enhanced 
backscattering, or from the dependence of the transmitted intensity on the length L of the 
disordered sample). For the quasi one-dimensional geometry of the present experiment, 
one expects V and I to obey the relation D = (1/3) VeI, with ve the energy transport 
velocity through the disordered medium. In the present case, the ratio 3D jl yielded a value 
ve = (5 ± 1) 10 7 m s _1 which is about an order of magnitude smaller than the phase 
velocity. This surprising result has been understood both for small concentration [0-0] and, 
more recently, also for strong concentration [ ]TU > |IT1 of the scatterers. (In the latter case, the 
resonant structure in the frequency dependence of T> disappears, and an overall decrease of 
V is observed 



|T2|1). Qualitatively speaking, the transport velocity is reduced because on its 
way through the medium, the energy is stored for some time in the Mie resonances. 

The present paper goes beyond previous theoretical work in two ways. First, we present 
a unified theoretical framework in which both average and fluctuation properties of the 
transmitted and reflected waves can be calculated on the same footing. This is in contrast 
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to approaches 0,|1O|,|TT[ which use the Bethe-Salpeter equation or a mean-field approxi- 
mation fTOlJTTH for the calculation of the transport velocity ve, and a diagrammatic impurity 
perturbation expansion for the intensity autocorrelation function. Second, we identify the 
energy dependence of the mean level density p(E) as the culprit for the observed deviation 
from standard behavior. A simple argument which yields the same result as the analytical 
derivation and illuminates the physical content may be helpful at this point. In the case of 
electrons [I], the Thouless formula g = E c p(E) connects the average conductance g with the 
Thouless energy E c = hV/L 2 and the mean level density. The presence of numerous reso- 
nances with equal resonance energies E\ leads to a local Breit-Wigner-shaped increase of 
p{E) near E\. Since g is not affected by the presence of the resonances, the Thouless formula 
implies that E c and, hence, the diffusion constant T> have a Breit-Wigner-like dip near E\. 
In Refs. |fL3| , |l4 | we have shown that within the framework of the supersymmetry formalism, 
the effective Lagrangeans for electrons and for classical light waves coincide. Hence, the 
argument just given applies likewise to scattering of light. With increasing concentration 
of scatterers, the dip in V widens and becomes less deep. Eventually, this results in an 
overall decrease of T> over a wide frequency interval. We emphasize that in the context of 
the supersymmetry approach, the relation V = (1/3)veI is not used explicitly. It is replaced 
by the Thouless relation. 

A summary of our results has been given in Ref. [|13|]. First results for electrons were 
published in Ref. 
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In the present paper we give a detailed derivation of our results, 
paying special attention to the necessary modifications of random matrix theory, and of the 
standard supersymmetry formalism. For classical waves, the latter is introduced in Ref. fll4 |. 

The random matrix model with resonances is introduced in section O. In section [TH], 
a particularly simple case is presented which paves the way for the full results given m 
section |V[ Section [V| contains the conclusions. We use Efetov's supersymmetric generating 
functional ]I3 in the version of ref. WJ\ . 



II. THE MODEL 

In Ref. |ll 5|| , we have given a detailed physical justification for the choice of the random 
matrix modeTused to describe diffusive scattering in the presence of resonances. The argu- 
ments presented there apply to the case of electrons. Meanwhile, we have shown that the 
non-linear sigma models for electrons and classical waves are identical, save for the source 
terms [13jT^]. For these reasons, we keep our presentation brief. In this section we always 



refer to electrons. It must be understood that all arguments apply equally to light (classical 

waves) 

In the absence of Mie resonances, a quasi one-dimensional disordered system of length 



L for diffusive scattering of electrons is modeled [HI as consisting of many slices of length 
lo <C L each. Within each slice, the Hamiltonian is modeled as a member of the GOE, 
the random-matrix ensemble with orthogonal symmetry. Neighboring slices are coupled by 
Gaussian-distributed uncorrelated random matrix elements. The strength of this coupling 
defines the diffusion constant. The first and the last slice are coupled to the channels, i.e., 
the asymptotic states in the leads. Elements of the scattering matrix connecting incident 
and outgoing channels in both leads define the transmission through the disordered region 
and, thus, the conductance. 

To account for the presence of Mie resonances, we retain the idea of a division of the 
quasi one-dimensional system into slices but modify the form of the Hamiltonian H. Within 
each slice, H is a matrix of dimension N + m where m is the number of Mie scatterers within 
each slice. This number is determined by the concentration of scatterers, and by the linear 
dimensions of the slice. Thus, H has the form 



H " ( V™ E,x I m + H„) (1) 

where T denotes the transpose. Here, -£/goe is a matrix of dimension N belonging to the 
Gaussian orthogonal ensemble (GOE) of random matrices, E\ x I m is an m-dimensional 



2 



diagonal matrix with elements E\ = huo\ corresponding to the presence of m Mie scatterers 
with equal resonance frequencies U\, I m is the m-dimensional unit matrix, H res is a random 
matrix of dimension m, also a member of the GOE which describes the coupling between res- 
onances, and V is a rectangular matrix which couples the m resonances to -£/goe- We assume 
the matrices V, H Tes , and H GOE to be statistically uncorrelated. The variance of the matrix 
elements of Hqoe is chosen as usual: The ensemble average of (Hcon) ^(Hgoe) ^v' is given by 
(\ 2 /N) (S fl y5 u y + 5^y5 v y). Here, A has the dimension of an energy and determines the av- 
erage level spacing ds (the same in each slice). For m = and with A (E) = \Jl — (E / (2\)) 2 , 

the dependence of ds on energy E is given by n\/(NA(E)) (semicircle law). At the end 
of the calculation, we take the limit iV — > oo. The ensemble average of (H Tes ) ^ u (H rcs ) ^ „/ is 
written analogously as (A 2 /m)(5^y5 u y + 5^y5 v y). The strength factor Ai is chosen in such 
a way that the interaction between resonances results in a lifting of the degeneracy which 
is of the order of the mean level spacing ds(0) so that Ai A. Too strong an interaction 
would wash out the resonance structure altogether. Finally, and without loss of generality, 
the rectangular matrix V can be taken to be diagonal. Indeed, V can be written in the 
form OiVdOJ where Vd is diagonal and rectangular, and where 0\ and O2 are orthogonal 
matrices of dimensions N and m, respectively. Transforming H with the orthogonal matrix 
O = diag(Oi,C>2), and using the orthogonal invariance of the GOE to absorb the matrices 
0\ into Hqoe and O2 into H res , respectively, we obtain a new Hamiltonian of the form 
of Eq. ([!]), but with V replaced by the diagonal rectangular matrix Vd- Since we assume 
all scatterers to be identical, we take all diagonal matrix elements of Vd to be equal. We 
assume that the diagonal matrix elements v scale with N in the same way as those of the 
matrices H GO e an d H Tes , so that v 2 = a\ 2 /N where a is independent of N. This completes 
the definition of the model. 

It is useful to relate Ai and v microscopically to the properties of the Mie resonance, 
and to the concentration /i of scatterers. By construction, the total coupling strength of 
each of the m resonant states is on average given by v 2 + A?. This quantity must, therefore, 
be proportional to the square of the matrix element for electromagnetic decay of the Mie 
resonance. Moreover, Ai must grow monotonically with increasing concentration of scatterers 
because an ever greater fraction of the light emitted from one resonance will be absorbed 
by another one, rather than scattered randomly. More precisely, consider a granule with a 
Mie resonance emitting a spherical wave. This wave hits an (arbitrary) second granule at 
distance s. For geometrical reasons, the cross section for absorption by excitation of the Mie 
resonance in this second granule is proportional to s~ 2 . With changing concentration, the 
distance s between any pair of granules changes as s 3 ~ m _1 . Hence, the coupling strength 
between any pair of granules (which we identify with \\/m) is proportional to m 2 / 3 . This 

yields Ai ~ m 5//6 . This consideration suggests that the coupling between granules should be 
dependent upon their distance, a fact not taken into account in Eq. (|l|) where all resonances 
are coupled in the same way. Thus, using a random band matrix instead of the GOE matrix 
H res may yield a better simulation of the actual situation. We shall find, however, that the 
entire effect of H res consists in widening the peak in the mean density of states. Therefore, 
it does not seem worthwhile to go through the trouble of dealing with a more complicated 
ensemble. On the other hand, this consideration shows why it is justified not to take into 
account the interaction of granules located in different slices. This is the approximation we 
will use. 



III. ONE SLICE WITH COUPLED RESONANCES 



For pedagogical reasons we focus attention in this section to the simplest case of a single 
slice described by the Hamiltonian H given by Eq. (|I|). In Ref. (TjJ, we have studied the 
case of a single slice with uncoupled resonances (H res = 0). Here, we generalize this study 
to the case of coupled resonances (H rcs 7^ 0). We calculate the mean level density and the 
level-level correlation function, using the supersymmetry technique of refs. fl6||rP 
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The advanced and retarded Green functions can be written as integrals over supervectors 



G%(E±u) = t*J d[*]^ex P [±i£(*)] 



(2) 



where the Lagrangean is given by 



£ = -*t ± w ± «j _ H) tf. 



The supervectors ^ are defined by 



^ = (S\S 2 , X , -**)■ 



(3) 



(4) 



The vectors S 1 and S 2 have N + m ordinary real components, whereas the vectors x have 
N + m anticommuting components. We define a source term J =diag in 
graded space, where j is a matrix of dimension (N + m) with entries J Ufl , and introduce the 
generating functional 



Z ± {E ±uj,J) = J d[^]exp 



±i(C - J#) 



(5) 



The Green functions are obtained as functional derivatives of Z with respect to J at J = 0, 



G 



1 5Z d 



A" 7 



4 5 J, 



(6) 



We use this expression to calculate the average of the product of a retarded and an advanced 
Green function taken at different energies, < G + G~ >. This quantity (the "two-point 
function" ) plays an important role in describing average properties of random systems, such 
as the level-level correlation function. Except for the dimensions of the supervectors ^ 
and of the a- and Q-matrices appearing below, and except for the occurrence of additional 
frequency variables, the 2fc-point function for any positive integer k is governed by an 
effective Lagrangean of the same type. 

The generating functional Z for the two-point function is given by 



Z(E,uj, J) = J d[^]exp 



-i^L 1/2 (E + ujL + ir)L — H — J) L l/2 m 



(7) 



We have introduced two 8(iV + m) -dimensional graded matrices, 
namely L =diag(J, J, - J, —J, J, J, -/, -/) and J =diag(j, j, -j, -j, j, j, -j, -j). Here, I 
denotes the iV + m-dimensional unit matrix. The quantity ^ now denotes a supervector 
with 8(N + m) components. All quantities are given in "advanced-retarded" notation (see 
Ref. flTI). We average over the two uncorrelated random ensembles H GOE and H res . The 
presence of these two ensembles causes the Hubbard-Stratonovich transformation to involve 
two supermatrices a and <j\. For the generating functional, we find 



Z(E,uj,J) = J d[er] exp 



N m 1 

■— trga 2 - — trgcr 2 - -Trglog(M — uL — J) 



where we define 



E-Xa -V D 
M= | E - Act 

-V D E-Ex- Ai<7i 




(8) 



(9) 
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Here, trg denotes the trace over an 8 x 8 matrix in graded space, and Trg the trace over an 
8(N + m) x 8(N + m) matrix in both graded and Hilbert space. The symbol d[er] stands for 
the differentials of all the independent variables in both a and o\. The dimensions of the 
block matrices in Hilbert space are indicated in Eq. (Q). After taking the limit N — > oo and 
keeping m fixed, variation with respect to a yields the standard saddle-point equation with 
a constant diagonal solution 



E 
2A 



±iA, 



(10) 



where A 
equation 



1 — (E/2X) 2 . Variation with respect to o"i produces a second saddle-point 

Ai 



0"! 



E — E\ — Xi&i — Tar' 



(11) 



where T = v 2 /X. Eq. ([TT|) defines the diagonal elements of o\ in terms of those of a. 
The saddle-point manifold is generated by the usual choice of signs in Eq. (|10"D , and by 
transforming the resulting diagonal matrices od and o\d both by the same matrix T. The 
latter is defined in Eq. (5.28) of Ref. ||17|| . The resulting matrices are denoted by Q and Qi, 
respectively. 

We integrate over the massive modes and expand the logarithm up to terms linear in uo 
and J, using the identity 



exp 



exp 



trg log 



.4 



BD~ 1 B T 





D 



(12) 



and the fact that traces over terms not containing either L or J vanish. Introducing the 
diagonal graded matrix A =diag(l, 1, 1, 1, —1, —1, —1, —1), we find 



Z(E,u,J) = J dfi(t) exp 



T trg[(g + A^ gi)(wA + J)] 



(13) 



The symbol d//(i) denotes the invariant measure, see Ref. fll7| . 

The expression in Eq. (|i~3D is formally equivalent to the standard one, obtained in the 
absence of resonances. This is seen by introducing a new diagonal supermatrix 



* mA 

and by defining a new Q-matrix Q* by 

Q* = Ta*T-\ 

This leads to 



Z{E,u,J) = J d/i(t) 



exp 



_^trg[(Q*(u;A + J)] 



(14) 



(15) 



(16) 



Because of the presence of the m resonances, the average density of states is not given 
directly in terms of the saddle-point solution <jd but is obtained by differentiation of the 
generating functional. This yields 



- — Imo 

A7T 



NA 



D 



\7T 



III 



XlTl 



Imcri . 



(17) 
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The integral over p must be equal to the dimension of the Hamiltonian, / p(E)dE = N + m. 
This can be verified by a straightforward calculation. To obtain a better understanding of 
the form of the density of states given by Eq. ([17]), we consider the two limits of very small 
and very large Ai. In the limit Ai — > oo (Ai S> T), the solution of Eq. ([11]) to zeroth order 
in T is given by 

E — E\ i , . 

<TlD = ~2\T + ' (18) 



where A' = yl — ((£? — £'i)/2Ai) 2 . For sufficiently large Ai, i.e. sufficiently strong coupling 

between the resonances, the coupling of the Hamiltonians Hqoe and H res can be neglected, 
and the density of states reduces to the sum of two semicircles 

NA mA' , s 

We see that the resonant structure is completely washed out. This was to be expected. 
In the opposite limit Ai — > (uncoupled resonances), it is easily seen that we retrieve our 
earlier result for uncoupled resonances fl5|l 

A(E)N TA/tt 



P{E) = ~J^~ + m (E- El - (E/2X)ry + (taJ > (20) 

where irX/ (NA) is the mean level spacing d$ in the absence of resonances. The second term 
in Eq. (^) has area m, peak height m/(7iTA), and Lorentzian shape. As expected, Eq. ( |20|) 
shows a resonance enhancement of the mean level density, centered at the energy E\ of the 
m resonances. The total width (rA) is determined by the strength v of the coupling to the 
random scatterers. We recall the relation v 2 = a\ 2 /N and find TA = aXA/N = (a/n)ds- 
This shows that T is of the order of the mean level spacing ds- For light scattering, we have 
a ^> 1. For formal reasons, we also consider the case a ~ 1 in this paper. 

The expression for the width w of the Lorentzian in Eq. ( P0|) can be derived from a 
different point of view with the help of the golden rule. With |r > a resonance state and 
\i > a state in the iV-dimensional subspace of Hilbert space comprising the GOE part of H, 
we find, averaging over the states labelled \i >, 

1 N 

w = 7r-J2<r\V\t> 2 6(E t - E x ) . (21) 

i 

Since X)f ' S(Ei — E\) = NA/ttX, we immediately find w = TA, i.e. the result we obtain 
using the supersymmetry formalism. 

Using the explicit form of the matrix Q* in Eq. (|15D, we can write the effective action in 
the exponent of Eq. ([16]) in the form 



z7r-^-trgT -1 ATA, (22) 

«eff 



where d e g(E) is the inverse of the average level density p(E) given in Eq. ([17]) . We note that, 
aside from the definition of d e g, the form of this expression is the same as in the absence 
of any resonances. Taking into account the fact that we can justify the second saddle-point 



equation only for lu <^ T, we conclude that (as in the case of uncoupled resonances JT5J) the 
only change in the level-level correlation function caused by the presence of resonances is 
due to the rescaling of the mean level density. Fig. 1 shows the result for a set of parameters 
given in the figure caption. 
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IV. CONDUCTANCE OF A QUASI ONE-DIMENSIONAL SAMPLE 



In this Section we consider a quasi one-dimensional system coupled to two perfect leads. 
As mentioned in Sec. O, we consider this system as consisting of K slices labelled % or j 
with i,j = 1,...,K. Within each slice, the Hamiltonian has the form given in Eq. (2). 
Neighboring slices (\j — i\ = 1) are coupled by random matrices Hjf v with Gaussian- 
distributed uncorrelated matrix elements with zero mean value and a second moment given 
by {p 2 /N){8^8 v y + 5^. u'Svy). For reasons given earlier, we assume that only the random 
parts of the Hamiltonians in neighboring slices are coupled. 

In the case of a single slice, we have seen that the mutual coupling of resonances leads 
only to a quicker broadening of the resonances than would be the case otherwise, and finally 
to the disappearance of the resonance enhancement effect. With the intention of simplifying 
the presentation we, therefore, now consider the situation where the resonances are not 
coupled even within one slice (Ai = 0). 

The coupling between the left (right) lead and the first (last) slice (j = 1 or K) is 
described by the real matrix elements W-[„ ||18|| . Here, the index a labels the channels in the 
leads. We assume that the dependence on energy of these matrix elements is negligible. We 
are interested in the influence of the resonances on the coupling to the channels. Therefore, 
we allow for a direct coupling between leads and resonances. Hence, the index [i runs from 
1 to N + m. Following Ref. Hl8| , we assume for j = 1,K the orthogonality relations 



N+m 

E Wiw£ = (JV + ro) {v{) 2 5 ab . (23) 

We use the many-channel Landauer formula 

= £(W| 2 + |S2T) (24) 

u.b 

which expresses the conductance g at zero temperature in terms of the elements S^ b R of the 
scattering matrix associated with the total Hamiltonian. Here R and L stand for the right 
and left leads, respectively. The scattering matrix has the form 

S c a t = 5 c % b - 2Z7T £ W^[D-% d u WL (25) 

where the inverse propagator D is a matrix of dimension K x (N + m) given by 

D% = E8^ -H% + ikY: W^W c av {5^cL + 8 iK 5 cR )5^, (26) 



where i,j = 1, 2, . . . , K and /j,, v = 1, 2, . . . , N + m. The indices (c, d) on D 1 in Eq. fl25|) 
stand for (1,1) if (c,d) = (L,L), for (1,K) if (c, <f) = (L,R), etc. The energy E equals the 
Fermi energy and eventually will be taken to have the value zero. 

To calculate the ensemble average of g, we follow Refs. |l7j and |18| and express the 



conductance in terms of a supersymmetric generating functional. We only indicate the steps 
where we deviate from the references just given. We define a generating functional Z for 
the conductance g, in a way similar to that presented in the previous section. 



Z(J)= Jvfi/]exp U^L 1/2 ( y E + iQ(j)+ir ] L-H_ + J( y j))L 1/2 y 



(27) 



The definitions of the source term J and of the matrix Q may be found in Ref. [|l^]. Suffice 
it to say that both J and Q are different from zero only for the first and the last slice (only 
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these two slices are coupled to the leads), and that Q essentially comprises the last term 
on the rhs of Eq. (p6|). After averaging and the Hubbard-Stratonovich transformation, the 
generating functional for the conductance g has the form 



Z(J) = J d[a] exp 



•^^trga^Mj) - ^IVglog(M(j) +iQ(j) + J(j)) 

i,3 3 



(28) 



As before, the indices i,j label the slices. The integration extends over K independent 
er-matrices. A matrix g = P^ 1 of dimension K (not to be confused with the conductance) 
is introduced, where P^ = (1/N)[5ij + (p 2 / X 2 )(5ij + i + 5 i+ ij)}. For each j, the matrix M(j) 
has the form given in Eq. @ with Ai = 0. 

Using the saddle-point approximation (valid for N ^> 1), we determine the mean level 
density and show that the resonant structure found for a single slice remains qualitatively 
unchanged. Omitting the source term, and taking the limit m/N — > 0, we obtain the 
saddle-point equation 



(Tit 



X N 



E - Xa(j) 



(29) 



We note that this equation is identical to the one which would obtain in the absence of all 
resonances. It was shown in Ref. [n| that the unique diagonal solution of Eq. (|29| ) has the 
form 



a D (j) = r j I 4 -iA j L. 



(30) 



Assuming weak coupling between the slices, we find for the total density of states in the 
entire sample in zeroth order in p 2 j X 2 



-tot 



X>(j), 

i=i 



(31) 



where p(j) is given by Eq. fl2"D| ) with r J substituted for E/2X, A- 7 for A. This result shows 
that the Lorentzian part in the average level density which is due to the contribution from 
the resonances, is preserved. The importance of this statement lies in the fact that the 
conclusions valid for a single slice apply likewise to the entire sample. 

To see this, we calculate the two-point function. We follow Ref. |18[ and parameterize 
a(j) as a(j) = (T^y l a D (j)T^ with 



T 3 



1 + ti2^21 



it{ 2 



-it j 



21 



1 + f t 3 
1 -|- t 2 l , '12 



(32) 



We emphasize that in spite of the presence of the m resonances in each slice, the form of 
the matrices a(j) is the same as in Ref. |18| (absence of resonances). We recall that the 
matrices Q and J differ from zero for the nrst and the last slice only. From the identity 
Eq. ([12] ) and the fact that Trglog(£' — Ei) = 0, we conclude that the logarithmic term in 
the exponent of Eq. (E3) vanishes unless j = 1 or j = K. Except for these two terms, 



the resulting form of the generating functional is identical to the one given in Ref. ||T8 
Moreover, the terms with j — 1 or K differ from zero only because of the presence of the 
matrices Q and J both of which involve the coupling to the channels. These results show 
that the presence of resonances does not affect the behavior of the average conductance in 
the bulk, but only the way it is coupled to the leads. It was shown in Ref. |L8| that for 
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sufficiently long samples (where the Thouless energy is smaller than the decay width defined 
by the coupling to the leads), the influence of the coupling to the leads is immaterial for 
the behavior of the average conductance. In practice, this condition is met for samples with 
lengths exceeding a few elastic mean free paths. Using these facts, we conclude that for 
typical samples, the average conductance is the same with and without resonances. Since 
the average conductance is calculated at fixed frequency, this conclusion holds likewise for 
electrons and for light (provided the frequency in the latter case is sufficiently high). This 
is also true for the calculation of the transmission coefficients given below. 

Before addressing the influence of the resonances on the source term and on the cou- 
pling to the channels, we turn attention briefly to the intensity autocorrelation function 

g(E)g(E + to). We do not describe the calculation of this quantity (a four-point function) 
in any detail. It is not difficult to see, however, that the effective action for this quantity 
also has the form given in Eq. fl28|), except for the following changes, (i) The dimensions 
of the matrices a(j) are doubled, (ii) For each value of j = 1, . . . , K, an extra term — ujA 
appears under the log arithm. Here, A is the sixteen-dimensional analogue of the matrix 
defined above Eq. (|T3|). 

What are the changes caused by the appearance of this extra term under the logarithm? 
We neglect the source term and the Q matrix. For each value of j, we expand the logarithm 
in powers of u, keeping only the linear term. (The term of zeroth order vanishes). After 
some simple algebra, we find that the result has a form analogous to expression (|22| ) , 

UJ K 

- in— ]T trgCP^AT^A . (33) 

"eff j = i 

It has to be borne in mind that the T's are now matrices of dimension sixteen. Here, d e g 
is defined as the inverse of the mean level density p per slice, see Eq. (^), but with r J 
substituted for E/2X, A J for A, see Eq. (|30|). We assume that p has the same value in each 
slice. Comparing our result with what would apply in the absence of resonances, we see 
that in each slice, 1/d is replaced by l/rf e fr- For Efetov's effective action (obtained in the 
continuum limit where the number of slices tends to infinity while their individual length 
shrinks to zero) this means that the frequency u is multiplied by the total level density 
p tot defined in Eq. ( |3Tj) . This is the claim made in Ref. ]TI| and in the Introduction to the 
present paper. In contradistinction to the average conductance, the four-point correlation 
function is influenced by the presence of resonances: The frequency u occurs only in the 
combination up tot . The result for T> is shown in Figure 1 of Ref. [|l3j. The argument can 
easily be extended to correlation functions of higher order. 

The conclusions drawn in the last three paragraphs contain the main results of this 
work. It remains to calculate the influence of the resonances on the coupling to the leads. 
For reasons given above, this may be considered a somewhat academic exercise. However, 
experiments on short samples would be able to test the predictions made in this way. 

The terms describing the coupling to the channels are part of the matrices Q and J under 
the logarithm in Eq. (p8|). We recall that both these terms are different from zero only for 
the first (j = 1) and last (j = K) slices. In expanding the logarithm, we have to keep terms 
of any order in iQ(j) but only of zeroth and first order in J. Here, we focus attention on the 
terms of zeroth order in J. (The calculation of the first-order terms proceeds analogously). 
After calculating the trace over /i in each term of the series by using the orthogonality 
relation Eq. ( p3|) for j = 1 (K) and resumming the result, we use the saddle-point solution 
for the first and the last slices, the parameterization of a given in Eq. (|32|) , and introduce 
the quantities T a = V a 2 /A. After some algebra, we find for j = 1, K 

-^Trglog(l + iT a 44), (34) 
where T a is the "sticking probability" in channel a, given by 
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(i | Wa A r 2 A r-f(g-gi) \ 2 / gWa r2 (g-gi)(i-^(f) 2 +#r y ' 

y A A (ij_£; 1 — ^.r) 2 +r 2 A 2 y y2A A A _|Lr) 2 +r 2 A 2 / 

The sticking probabilities T a are related to the matrix elements of the average scattering 
matrix by T a = 1 — |S^| 2 . This means that T a measures that part of the flux which is not 
re-emitted instantaneously into the same channel but which experiences multiple scattering 
in the disordered sample. Eq. (|35|) very clearly shows that the sticking probabilities are 
affected by a direct coupling between resonances and channels. 



V. CONCLUSIONS 



It was shown in Ref. [|T8j that the conductance of a short sample is proportional to the 
product of two terms, each term being given by the sum of the sticking probabilities of the 
channels in one of the two leads. Thus, the conductance of a short sample is determined by 
the coupling to the leads, whereas the conductance of a long sample is entirely determined 
by the matrix appearing in Eq. This implies that as long as there is no direct 

coupling between resonances and leads, the conductances of both short and long samples 
are unchanged. When the resonances are coupled to the channels, however, the conductance 
of a short sample can be changed drastically at frequencies near the resonance frequency, 
whereas the corresponding change of the conductance of a long sample is expected to be 
small. It would be interesting to search for this effect experimentally. 

On the other hand, the intensity autocorrelation function is always modified by the 
presence of resonances, both for short and for long samples. This fact follows directly from 
our calculation of the diffusion constant T> (which shows a resonance dip) but is most easily 
seen by using an argument based on the Thouless relation. This argument was given in the 
Introduction and is not repeated here. The argument is fully corroborated by our analytical 
results. 

It was emphasized above that the effective Lagrangean for electrons and for scalar waves 
(light) is the same, both in the absence and in the presence of Mie resonances, provided that 
the usual condition kl 3> 1 is met. Here, k is the wave number and I the elastic mean free 
path. In view of remarks in the literature that the strong effect of Mie resonances upon light 
scattering is absent in the case of electrons, we wish to clarify this statement. In the case of 
scalar waves (light), there are two different and unrelated effects which contribute towards 
the energy dependence of the potential. The first of these is due to the dependence on k of the 
term k 2 e(r) in the scalar wave equation. Here, e is the space-dependent dielectric constant. 
The disorder potential V in the Schrodinger equation has no analogous dependence on wave 
number or energy. However, the change with k of k 2 e(r) is small on the scale of the mean 
level spacing. This is why the non-linear sigma models for electrons and for scalar waves are 
identical, even though the Ward identities for the two wave equations differ. Second, the 
presence of Mie resonances causes the drastic modification of the intensity autocorrelation 
function for light scattering discussed above. The equality of the nonlinear sigma models 
for scalar and for Schrodinger waves implies that a correspondingly large modification of the 
intensity autocorrelation function would also occur for electrons if it were possible to design 
a system in which electrons are scattered both diffusively and via Mie-like resonances. To 
the best of our knowledge, such a system does not exist at present, and the stated analogy 
makes a prediction which cannot be checked experimentally. 

To compare with experiment, we note: (i) For classical waves in a uniform medium in d 

dimensions, the density of states is proportional to e d l 2 where e is the dielectric constant of 
the medium. Increasing the number of resonating spheres, we increase the effective dielectric 
constant of the medium. Hence, A should decrease with increasing /x. (ii) As shown in Sec. [TI|, 
Ai is expected to increase monotonically with /x. Both these points affect the dependence 
on concentration of d e g and, therefore, of V ~ d cS . With increasing /x, point (i) causes a 
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decrease of T> far away from the resonance dip, and point 
Fig. 1 of Ref. p| 



ii) a widening the resonance dip. 
shows T> versus energy for different concentrations. In agreement with 
the experimental results ||12], we find a deep dip in the diffusion constant at the resonance 
energy. It is caused by the second term in Eq. (^) and is particularly pronounced at low 
concentrations of the resonant scatterers. As the concentration increases, the dip is smeared 
out. At the same time, the value of V outside the resonance dip decreases. 
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FIGURES 




FIG. 1. The exact normalized level-level correlation function for m = 5, V = 5d, E 
and rescaled correlation function for a pure GUE given by sin 2 : (2ttuj / 'd e f /) / '(2-itlj / 'd e f /) 2 . 



13 



